Hydrodynamic theory of surface excitations of three-dimensional topological insulators 
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Edge excitations of a fractional quantum Hall system can be derived as surface excitations of an 
incompressible quantum droplet using one dimensional chiral bosonization. Here we show that an 
analogous approach can be developed to characterize surface states of three-dimensional time reversal 
invariant topological insulators. The key ingredient of our theory is the Luther's multidimensional 
bosonization construction. 
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Topological insulators are materials with insulating 
bulk and topologically protected metallic edge or surface 
states*^— Two-dimensional TIs are also called as quan- 
tum spin Hall systems. Edge states of QSH systems are 
analogous to quantum Hall edge states. However there 
is an important distinction^: edge states of the QHE 
with up and down spins propagate in the same direction, 
whereas QSH edge states with opposite directions of spin 
counter-propagate. This is the consequence of the time 
reversal invariance of QSHE, which is broken in conven- 
tional QH systems. Therefore the former are called chiral 
and the latter are called helical edge states. 

There are two classes of non-trivial time reversal in- 
variant TIs in three-dimensions which are called strong 
and weaker— While weak TIs are layered 2D QSH states 
(in the sense that these two systems can be adiabati- 
cally coupled to each other) , strong TIs are purely three- 
dimensional. Surface states of these TIs are massless 2D 
Dirac fermions. Weak and strong TIs are distinguished 
by the number of Dirac cones on the surface: strong TIs 
have an odd number and weak TIs have an even number 
of Dirac cones on the surface. Gapless modes of strong 
TIs are robust and insensitive to weak interactions and 
disorder. 

Surface states of three-dimensional time reversal in- 
variant TIs are spin filtered, such that (s(— k)) = 
— (s(k)), which means that spin density and charge cur- 
rent are coupled (see, e.g., Ref.[l]). This suggests that 
three dimensional strong TIs realize quantum spin Hall 
effect in every radial direction. "Radial direction" here 
means that one considers tomographic projection of the 
surface states on certain direction. This tomographically 
projected state in certain cases can be viewed as a two 
dimensional quantum spin Hall system. In this paper we 
explicitly show how to construct such mapping mathe- 
matically. In constructing this mapping we assume that 
a hydrodynamic theory of the QSHE edges can be de- 
veloped in the similar way as for QHE edges^r— Sev- 
eral interesting physical quantities characterizing three- 
dimensional TIs are introduced. 

It should be made clear that many of the ideas pre- 
sented here taken separately are not new. We have only 
combined them in a single picture. We also note that 
this paper is illustrative rather than rigorous. 

A good description of the hydrodynamic theory of the 



FQHE edge states can be found in [12j. Here we briefly 
mention the main points. Then we modify this theory to 
describe surface states of strong topological insulators. 

Suppose that the FQHE states are incompressible ir- 
rotational liquid without bulk excitations. Then the only 
low lying excitations are the surface deformations of the 
quantum droplet. The droplet is confined by a smooth 
potential well. The electric field of the potential well will 
generate a persistent current along the edge fluctuating 
part of which is given by 
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where v is the filling fraction, e z is the unit vector along 
the z axis, h(x) is the displacement of the edge from its 
equilibrium value, x is the coordinate along the edge. 

One-dimensional density of the edge wave p{x) is re- 
lated to the displacement of the edge h(x) through p(x) = 
nh(x), where n = v^- is the two dimensional electron 
density in the bulk. Then continuity equation reads 



d t h — vd x h = 



(1) 



where v = cj|. This means that the electrons at the edge 
drift with the velocity v. 

Hamiltonian (energy) of the edge waves is given by 



H 



hpEdx 



7TD 
V 



p dx 



(2) 



It is easy to quantize this Hamiltonian. Rewriting ([T]) 
and ^ in the momentum representation and identify- 
ing pk with the coordinate variable one finds that the 
corresponding canonical momentum is given by pk = 
2nip^k/vk (zero mode with k = is excluded from the 
theory due to incompressibility of the liquid) . From the 
commutation relations \pk,Pk'] = *<5fcfc' one obtains the 
Kac-Moody algebra for the currents 
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This theory provides complete description of low lying 
excitations of the Laughlin stated 

Now we want to employ this theory to QSH systems. 
The easiest way to understand QSHE is to consider the 
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case when spin s z is conserved Then one can define two 
sectors with spin up and down. Such decomposition is 
possible even when spin is not conserved^ Then one can 
define the Chern numbers n^j, for the spin up and down 
sectors in the usual manneri 15 ' 16 In a time reversal invari- 
ant system the total Chern number + nj, = is zero. 
However the difference of the two Chern numbers in gen- 
eral is not zero and one can define the spin Chern number 
according to C s = \{n^ — 714,). Loosely speaking, C s de- 
termines the number of gaplcss edge modes in the system. 
It was shown that the edges of a QSH system with even 
C s are localized by disorder, while edges of a QSH sys- 
tem with odd C s are robust against small time reversal 
invariant perturbationsi 17 i 18 The case n-f- = — = 1 is 
the simplest and to describe edge states of such a system 
one needs two uncoupled incompressible liquids p\ and 
P2 with the edge state Hamiltonians @, which are re- 
lated by time reversal symmetry. Edge excitations travel 
in opposite directions with the velocity u = eE/hn. Here 
the index 1 refers to spin up and the index 2 refers to spin 
down sectors. 

In analogy with the above case, we assume that the 
low lying excitations of a three-dimensional topological 
insulator are the surface excitations of two uncoupled 
incompressible quantum liquids confined by a potential 
well, with the following Hamiltonian: 

H = j ^p eE[h 2 1 (x,y) + h 2 2 {x,y)]dxdy (4) 

where po is the density of the electronic liquid, x = 
x and y are coordinates along the surface. Two- 
dimensional density of the surface waves is pi(x,y) = 
pohi(x, y), i = 1,2. Incompressibility of the liquid means 
that J hi(x,y)dxdy = 0. It will turn out that after mak- 
ing some assumptions about the dynamics of the fields 
hig, this case corresponds to a strong TI with a single 
Dirac cone on the surface. 

The essence of Luther's approach is to consider 
a D-dimensional space as a set of one-dimensional 
"rays"ii*^— Following these ideas, we write the surface 
state Hamiltonian (j4|), as a sum over tomographic pro- 
jections 

h= L de I \ neE + dc (5) 

integrating over the right hemi-circlc 1Z = {—it/ 2 < 9 < 
7r/2}. Here the surface displacement of the tomographi- 
cally projected system hg(^) is defined as 

M?) = J °° (0^) ' dk J co ^-OH^,y')dx'd y ' 

k = \k\k,k = (cos 9, sin 9), £ = k ■ x, £' = k • x'; since 
z coordinate remains intact, we can consider the planes 
(£, z) which are labeled by 9. We call these planes as to- 
mographic planes or tomographic projections of the ini- 
tial three-dimensional system, n is the two-dimensional 



density of the electronic liquid in the tomographic plane, 
which is determined from the consistency of the two de- 
scriptions. Electric field of the confining potential well 
is the same in both cases^ It is easily verified that 
tomographic projections are also incompressible liquids. 
Now we make the following assumption that these tomo- 
graphic planes are QSH systems. From the topological 
band theory we can be sure that this is correct at least 
for three values of the parameter 9 (see the discussion at 
the end of the paper) . The fact that this is correct for all 
8 will be justified below, because this is the only assump- 
tion that leads to the desired result: single Dirac cone 
on the surface which will be obtained after fermioniza- 
tion of the model. Then one finds that excitations have 
linear spectrum luu = u\k\ 1 where u = eE/hn. Thus one 
can relate the unknown parameter of the theory n to the 
parameters of the three-dimensional theory. 

The choice of the range of 9 is not unique and this 
reflects the fact that splitting of the full Hilbert space 
induced by the time reversal operation is not unique. 
When time reversal symmetry is preserved such split- 
ting is necessary to obtain a non-trivial Chern number 
of a 2D system, because the Chern number of the whole 
Hilbert space vanishes (in 3D one needs to consider cer- 
tain 2D sections of the Brillouin zone and further split 
them using time reversal operation). 

One can define Fourier components of hig(£): 

ke(k) = / M£)e~ ife ^ (7) 

Since by our assumption hig(k) are the edge modes 
of a QSH system, they must satisfy the equations of 
motion (continuity equations) dthig{k) — iukh\g(k) 
and dth,2e(k) = —iukh2e(k). From © one also has 

h w (±\k\) = v/^/ii(±k). This finally leads to 

<9 t /ii(±k) = ±iu\k\h 1 (±k) (8) 
d t h 2 (±k) = Tm\k\h 2 {±k) (9) 

when k is in the right hemi-circle. In the second- 
quantized form the Hamiltonian Q is 

H B = J2u\k\(aia k + pt0 k ) (10) 

where 

e.g, ak is related to hi(k) through 

«k = V ' po/n\k\h-i{-k), al = ^/p /n\k\hi(k) when 
k G 72.. The Hamiltonian (fT0"|) is half of the massless 
Klein-Gordon model, exactly what is needed to construct 
massless two-component Dirac fermion^ This is similar 
to that chiral bosons in one-dimension are the half of the 
ordinary bosons. One can associate the following bosonic 
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fields with this model 



S 
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' 2^J /„ e ~ ak/2dk 



x (a k e ife(fe - x) + /i.e.) (11) 
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e~ ak/2 dk 
x (^e ifc ^ x ) + /i.e.) (12) 



a is the cutoff which should be taken to zero at the end; 
S is the surface area. These fields do not correspond 
to any local observables. However, Luther showed that 
appropriate functions of these fields do. 

Here we review some of the details of the Luther's 
construction for completeness £^ Suppose we have a 
fermionic Hamiltonian 



Hp = u a\ik ■ (?)ak 



(13) 



This is a Hamiltonian of masslcss Dirac fcrmions. It can 
be diagonalizcd by a transformation 



cra k 



(14) 



where V = (— sin#, cos 9). Fermi operators transform 
according to 



a„ = e- MVff / 4 a k 



(15) 

The diagonal Hamiltonian is 

H' F = }^u\k\a'^a z a{ s . (16) 

k 

The boson representation is given by 

V>(M-x) = e- iwV ' ff/ V(M • *) (17) 



(Klein factors necessary to ensure anticommutation rela- 
tions arc omitted for simplicity; for details see 19H20| ). 
The usual Fermi fields arc given by 



^(x) 



Ti 



d9ijj(9, k ■ x) 



(19) 



The representation (|IT]) . ([T2"j) . (|17lll9|) is constructed in 
such a way that the correlation functions of free fermionic 
fields are correctly reproduced. Note a crucial point that 
in the definition ([Ti?]) the integration is only over half of 
the whole circle. Thus we see that the Hamiltonian ((4]) 
together with the equations of motions ([5]) and © (or 



alternatively the Hamiltonian where /iie(£) are edge 
states of a QSH insulator) is equivalent to masslcss Dirac 
Hamiltonian. This corresponds to strong TI with a sin- 
gle Dirac cone on the surface. These calculations confirm 
our initial intuition. 

There is only one Z^ invariant in two-dimensions^ It 
can be defined as C s mod 2<2^ In three-dimensions there 
are four Z2 invariants Three of them are equivalent to 
invariants of two-dimensional topological insulators and 
are defined as invariants of some sections of the Brillouin 
zone. The fourth topological invariant i>q is purely three- 
dimensional. TIs with ^o = are called weak and TIs 
with vq = 1 are called strong. This invariant determines 
the number of Kramers degenerate Dirac points enclosed 
by the Fermi surface. 

In any time reversal invariant system with spin orbit 
interactions there are two-dimensional Dirac points in 
the surface spectrum. Therefore the above considerations 
should be clarified. Single Dirac fcrmion on the surface 
already means that the insulator is a strong topologi- 
cal insulator. We will show this directly and in parallel 
discuss the relation of our picture to the conventional 
theory of topological insulators. To make connection of 
this picture with the band topology, we use the simple 
argument for counting the topological invariants due to 
Roy2&, which is quoted below. For us the important as- 
pect of this work is how the Z2 invariants of certain planes 
in the Brillouin zone, such as p x — p y , can be computed. 

Let represent the Brillouin zone by a cube {— tt < 
Px,PyiPz < tt} and the Zi invariants associated with the 
planes p x = 0, p x = tt, p y = 0, p y = tt by v\,v\,v-i,Vi 
respectively. Then the Z^ invariant of the plane p x = p y 
equals v\ + vi and the Z^ invariant distinguishing strong 
topological insulators from weak topological insulators 
equals v§ = v\ — v\ = V2 — vi- Any 3D topological insu- 
lator with time reversal invariance can be characterized 
by four invariants, which may be chosen to be v\, ^2, ^3 
and vq- 

Above we assumed that every tomographic plane sup- 
ports QSHE. In fact, it is sufficient to consider only three 
planes (then for the rest this would be satisfied automat- 
ically). Let these planes be p x = n, p y = tt, p x = p y . We 
will consider for concreteness the Dirac point (tt, tt, 0) 
and the surface states on the (a;, y) plane having small 
2D momentum k around p x = ir,p y = tt (direction of k 
is given by the angle 9). In this case 9 = corresponds 
to the plane p y — tt, 9 = tt/2 to the plane p x = tt, and 
9 = tt /A to the plane p x = p y . Suppose that Z2 invariant 
of each of this planes is odd, thus corresponding to non- 
trivial insulator. Then we have v\ = = v\ + i>2 = 1, 
which gives v = 1. This corresponds to strong topolog- 
ical insulating phase. If there is no such a point in the 
Brillouin zone, for which all three Zi invariants are 1, 
then the insulator is in the weak topological insulating 
phase with = 0. (see also Fig. 3 in the Rcf. 27]) 

In general, it isn't possible to define a Z2 invariant of a 
plane with arbitrary 9 using topological band theory. It 
seems that this is possible only for planes such as tan 9 — 
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m/n, where m and n are two co-prime integers. In the 
case we have considered, all such invariants are equal to 
1. 

Recently, bosonization approach was applied to topo- 
logical insulators also in the work (28|, but in a differ- 
ent context. In passing we also note that surface excita- 
tions of certain 3D topological superconductors^ can be 
viewed, in every radial direction, as edge states of a 2D 
topological superconductor in the same class. It would 
be interesting to explore this case as well. 

In summary, we have shown that low lying excitations 
of a strong TI with a single Dirac cone on the surface 
are the surface deformations of a droplet of incompress- 
ible quantum liquid. These excitations have very unusual 
form, however they have a simple meaning when one con- 



siders tomographic projections of this liquid (as defined 
in the text): they are two chiral waves propagating in 
opposite directions. Thus, surface excitations of a strong 
TI with a single Dirac cone on the surface are the sum 
of QSH edge states. However, this is not true for the 
entire topological insulator, i.e., a strong TI can not be 
presented as a sum of 2D TIs. This paper can be consid- 
ered as an another physical illustration of the fact, that 
3D topological insulators can be characterized by 2D in- 
variants. 
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versions of this manuscript and Yu.E. Lozovik for draw- 
ing my attention to topological insulators. 
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Although Luther considered three-dimensional fermions, it 
is easy to extend his results to two dimensions. 



